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Abstract 

Beauty photoproduction in dijet events has been measured at HERA with the 
ZEUS detector using an integrated luminosity of 126 pb _1 . Beauty was identified 
in events with a muon in the final state by using the transverse momentum of 
the muon relative to the closest jet. Lifetime information from the silicon vertex 
detector was also used; the impact parameter of the muon with respect to the 
primary vertex was exploited to discriminate between signal and background. 
Cross sections for beauty production as a function of the muon and the jet 
variables as well as dijet correlations are compared to QCD predictions and to 
previous measurements. The data are well described by predictions from next- 
to-leading-order QCD. 
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1 Introduction 



The production of beauty quarks in ep collisions should be accurately calculable in pertur- 
bative Quantum Chromodynamics (QCD) since the large mass of the b quark provides a 
hard scale. Therefore it is interesting to compare such predictions to results using photo- 
production events where a low-virtuality photon, emitted by the incoming lepton, collides 
with a parton from the incoming proton. Previous photoproduction analyses presented 
by ZEUS used semi-leptonic decays into muons [1,2] and electrons [3,4], and found agree- 
ment with next-to-leading-order (NLO) QCD calculations and predictions based on the 
kr factorisation approach [5]. Leptons from beauty decays were distinguished from other 
decay leptons and background by their large transverse momentum, p™ 1 , relative to the 
axis of the jet with which they are associated. 

For the analysis presented here, events with at least two jets (jj) were selected and b 
quarks were identified through their decay into muons with large p™ 1 and large impact 
parameter, 5, defined as the distance of closest approach of the muon with respect to 
the beam position. The latter was facilitated by the ZEUS silicon microvertex detector 
(MVD) [6]. The impact parameter is large for muons from b decays due to the long 
lifetime of B hadrons. A combination of the p™ 1 and 5 methods was also used by the 
HI collaboration [7] and good agreement was found with both the ZEUS results and the 
NLO QCD prediction, except for an excess at low px of the muon (p^) and jet (p 3 T ). The 
measurement presented here covers a kinematic region extending to lower pj, than the 
previous ZEUS and HI jet measurements [1,7]. The cross section for beauty production 
has also been measured in pp collisions at the SppS [8] and Tevatron colliders [9] and in 
77 interactions at LEP [10,11]. Most results are in good agreement with QCD predictions 
but large discrepancies are observed in some [10] of the results from 77 interactions at 
LEP. 

The dijet sample of beauty photoproduction events was also used to study higher-order 
QCD topologies. At leading order (LO), the two jets in the event are produced back- 
to-back in azimuthal angle, such that A0- 7 - 7 = jl — j2 = tc. Additional soft radiation 
causes small azimuthal decorrelations, whilst A^- 7 significantly lower than ir is evidence 
of additional hard radiation. Dijet correlations have been previously measured at ZEUS 
in inclusive-jet and charm photoproduction at high transverse energies [12-14]; the con- 
clusions from both are the same. Deviations from the NLO QCD predictions were found, 
especially in regions which are expected to be particularly sensitive to higher-order ef- 
fects, i.e. at low A0-". In this paper, the cross section versus A0-" is presented for beauty 
photoproduction. These and other cross sections are compared to NLO QCD predictions 
and Monte Carlo models. 
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2 Experimental set-up 



The analysis was performed with data taken in 2005 when HERA collided electrons with 
energy E e = 27.5 GeV with protons of energy E p = 920 GeV, resulting in a centre-of-mass 
energy of \/s = 318 GeV. The results are based on an e~p sample corresponding to an 
integrated luminosity of 125.9 ± 3.3 pb _1 . 

A detailed description of the ZEUS detector can be found elsewhere [15]. A brief outline 
of the components that are most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central track- 
ing detector (CTD) [16] and the MVD [6]. These components operated in a magnetic field 
of 1.43 T provided by a thin superconducting solenoid. The CTD consisted of 72 cylin- 
drical drift chamber layers, organised in nine superlayers covering the polar-angl^] region 
15° < 6 < 164°. The MVD consisted of a barrel (BMVD) and a forward (FMVD) section 
with three cylindrical layers and four vertical planes of single-sided silicon strip sensors 
in the BMVD and FMVD respectively. The BMVD provided polar-angle coverage for 
tracks with three measurements from 30° to 150°. The FMVD extended the polar-angle 
coverage in the forward region to 7°. After alignment the single-hit resolution of the 
BMVD was 25 fim and the impact-parameter resolution of the CTD-BMVD system for 
high-momentum tracks was 100 /xm. 

The high-resolution uranium-scintillator calorimeter (CAL) [17] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a(E)/E = 0.18/a/E for electrons and 
a(E)/E = 0.35/ VE for hadrons (E in GeV). 

The muon system consisted of rear, barrel (R/BMUON) and forward (FMUON) tracking 
detectors. The B/RMUON consisted of limited-streamer (LS) tube chambers placed 
behind the BCAL (RCAL), inside and outside a magnetised iron yoke surrounding the 
CAL. The barrel and rear muon chambers cover polar angles from 34° to 135° and from 
135° to 171°, respectively. The muon system exploited the magnetic field of the iron yoke 
and, in the forward direction, of two iron toroids magnetised to ~ 1.6 T to provide an 
independent measurement of the muon momentum. 

The luminosity was measured using the Bethe-Heitler reaction ep —> ejp by a luminosity 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 
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detector which consisted of independent lead-scintillator calorimeter [18] and magnetic 
spectrometer [19] systems. The fractional systematic uncertainty on the measured lumi- 
nosity was 2.6%. 

3 Data selection 

The data were preselected by the ZEUS online trigger system [15,20,21] to contain two 
high-^T jets and/or a muon candidate. 

The hadronic system (including the muon) was reconstructed offline from energy-flow ob- 
jects (EFOs) [22] which combine the information from calorimetry and tracking and which 
were corrected for dead material and for the presence of muons. Jets were reconstructed 
offline from EFOs using the kx algorithm [23] in the longitudinally invariant mode [24]. 
The i^-recombination scheme, which produces massive jets whose four-momenta are the 
sum of the four-momenta of the clustered objects, was used. 

Muons were reconstructed by matching a track found in the CTD and the MVD with a 
track found in the inner chambers of the B/RMUON. Muons were associated with jets 
using the kx algorithm; if the EFO corresponding to a reconstructed muon was included 
in a jet then the muon was considered to be associated with the jet, which will from now 
on be referred to as the muon-jet. 

Events with one muon and two jets were selected by requiring: 

• > 1 muon with pseudorapidity —1.6 < if 1 < 1.3, and transverse momentum pj, > 
2.5 GeV (this cut was lowered to p£ > 1.5 GeV for the measurement of the differential 
cross section with respect to pj) ; the muon track was required to have at least 4 MVD 
hits; 

• > 2 jets with pseudorapidity [^l < 2.5, and transverse momentum p° T > 7 GeV for 
the highest-p^ jet and p° T > 6 GeV for the second-highest-p^ jet; 

• that the muon was associated with a jet with p 3 T > 6 GeV which is not necessarily 
one of the two highest-pr jets. To ensure a reliable p™ 1 measurement (see Section HJ), 
the residual jet transverse momentum, calculated excluding the associated muon, was 
required to be greater than 2 GeV; 

• no scattered-electron candidate [25] found in the CAL; 

• 0.2 < ?/jb < 0.8, where t/jb is the estimator of the inelasticity, y, measured from the 
EFOs according to the Jacquet-Blondel method [26]. 

The last two cuts suppress the contributions from neutral current deep inelastic scattering 
events and from non-ep interactions. The total efficiency of all these selection cuts was 
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~23%. A sample of 7351 events remained for pt^ > 2.5 GeV and 14172 events remained 
for p£ > 1.5 GeV. 



4 Signal extraction 

To evaluate detector acceptance and to provide the signal and background distributions, 
Monte Carlo (MC) samples of beauty, charm and light-flavour (LF) events were gen- 
erated using Pythia 6.2 [27-29], corresponding respectively to 9, 4.5 and 1 times the 
luminosity of the data. The generated events were passed through a full simulation of the 
ZEUS detector based on Geant 3.21 [30]. They were then subjected to the same trigger 
requirements and processed by the same reconstruction programs as the data. 

Due to the large 6-quark mass, muons from semi-leptonic beauty decays tend to be pro- 
duced with high transverse momentum with respect to the direction of the jet containing 
the B hadron. The pf. el variable can therefore be exploited to extract the beauty signal; 
it is defined as: 



\p* x (pi -p**) 



\p> - pA*| 

where p^ is the muon and p J the jet momentum vector. An underestimation of the tails 
of the p™ 1 distribution in the background MC was corrected as described in a previous 
publication [1]. The correction is applied in bins of p™ 1 to the LF MC sample. Half of 
the correction was also applied to the charm MC sample. The p™ 1 distribution for the 
selected data sample is compared to the MC simulation in Fig. QJa). 

Muons from semi-leptonic beauty decays tend to be produced at a secondary vertex, 
displaced from the primary vertex, because of the long lifetime of B hadrons. The signed 
impact parameter 5 is calculated with respect to the beam position in the transverse 
plane (beam-spot). The beam-spot position was calculated every 2000 events as described 
elsewhere [31]. The sign of 5 is positive if the muon intercepts the axis of the associated 
jet within the jet hemisphere; otherwise 6 is negative. 

The impact-parameter resolution in the MC was corrected [32], simultaneously taking 
into account the residual effects of multiple scattering and of the tracking resolution. 
The correction was extracted from inclusive-jet data by fitting the impact-parameter 
distribution with a double convolution of a Gaussian and a Breit-Wigner function. The 
widths of these functions were tuned taking into account the pt dependence of the size of 
the correction. 

Figure [H(b) shows the distribution of the reconstructed muon 5 compared to predictions 
from the Pythia MC model for beauty, charm and LF events which were corrected as 
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described above. The 5 distribution for LF events, which is symmetric and peaked at zero, 
has a finite width which reflects the impact-parameter resolution. Whereas for beauty 
events, and to a lesser extent for charm events, the 6 distribution is asymmetric towards 
positive 5. 

The fractions of beauty (a bb ) and charm (a cE ) events in the sample were obtained from a 
three-component fit, f^, to the measured two-dimensional distributions of p™ 1 and 5: 

U = G rf + <W? + (! - a bl ~ ^)/m LF > ( 2 ) 

where f bb , f^ c and f^f are the MC predicted shapes for beauty, charm and light flavour 
events. The fit used the minimum-x 2 method and included MC statistical uncertainties. 
For differential cross sections the fit was repeated for each bin. 

As an illustration, Fig. [2] shows the 68% probability contours from the two-dimensional 
fit described above and from one- dimensional fits carried out using p£ cl or 5 alone for the 
data sample with p£ > 2.5 GeV. Only statistical errors were taken into account. The two 
variables give complementary information. The p™ 1 fit alone is able to distinguish the b 
component from charm and LF but not to separate these two background components, 
while the 6 fit gives a good determination of the total heavy quark fraction. The 8 fit 
also provides a strong anti-correlation between the fractions of beauty and charnjfl In 
the previous analysis [1], which used the p™ 1 method alone, the charm contribution was 
constrained to the charm cross section obtained from other measurements. This is not 
necessary here. 



5 Theoretical predictions and uncertainties 

The measured cross sections are compared to NLO QCD predictions based on the FMNR 
[33] program. The parton distribution functions used for the nominal prediction were 
GRVG-HO [34] for the photon and CTEQ5M [35] for the proton. The 6-quark mass was 
set to rrib = 4.75 GeV, and the renormalisation and factorisation scales to the transverse 

mass, /i r = /i/= rriT= ^ \ ((Pt) 2 + G°t) 2 ) + m h where is the transverse momentum 

of the b (b) quark in the laboratory frame. Jets were reconstructed by running the kr 
algorithm on the four-momenta of the b and b quarks and of the third light parton (if 
present) generated by the program. The fragmentation of the b quark into a B hadron 
was simulated by rescaling the quark three- momentum (in the frame in which p b z = —p b z , 

2 Due to the correlations between pj? 1 and 6 the combined contour is not completely contained within 
the overlap of the two individual contours. 
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obtained with a boost along Z) according to the Peterson [36] fragmentation function 
with e = 0.0035. The muon momentum was generated isotropically in the B-hadron 
rest frame from the decay spectrum given by Pythia which is in good agreement with 
measurements made at B factories [37]. 

The NLO cross sections, calculated for jets made of partons, were corrected for jet hadroni- 
sation effects to allow a direct comparison with the measured hadron-level cross sections. 
The correction factors, Chad, were derived from the MC simulation as the ratio of the 
hadron-level to the parton-level MC cross section, where the parton level is defined as 
being the result of the parton-showering stage of the simulation. 

To evaluate the uncertainty on the NLO calculations, the 6-quark mass and the renormal- 
isation and factorisation scales were varied simultaneously to maximise the change, from 
m& = 4.5 GeV and /i r — (J,f= rriT/2 to m& = 5.0 GeV and fi r = fif = 2mr, producing a 
variation in the cross section from +34% to —22%. The effect on the cross section of a 
variation of the Peterson parameter e and of a change of the fragmentation function from 
the Peterson to the Kartvelishvili parameterisation was found in a previous publication [1] 
to be of the order of 3%. The effects of using different sets of parton densities and of 
a variation of the strong coupling constant were found to be within ±4%. These effects 
are negligible with respect to that of a variation of the fe-quark mass and the renormal- 
isation and factorisation scales and are therefore not included. The uncertainty due to 
the hadronisation correction was also found to be negligible with respect to the dominant 
uncertainty. 

The measured cross sections are also compared to the predictions of the Pythia 6.2 MC 
model scaled to the data. The predictions of Pythia were obtained [28] by mixing flavour- 
creation processes (jg —> bb, gg — + bb, qq — > bo) calculated using massive matrix elements 
and the flavour-excitation (FE) processes (bg —>■ bg, bq —>■ bq), in which a heavy quark is 
extracted from the photon or proton parton density. The FE processes contribute about 
27% to the total bb cross section. The small (~ 5%) [1] contribution from final-state gluon 
splitting in parton showers (g — > bb) was not included. The parton density CTEQ5L [35] 
was used for the proton and GRVG-LO [34] for the photon; the 6-quark mass was set to 
4.75 GeV and the fe-quark string fragmentation was performed according to the Peterson 
function with e = 0.0041 [36]. 

6 Systematic uncertainties 

The main experimental uncertainties are calculated as follows (the resulting uncertainty 
on the total cross section is given in parentheses) [38]: 
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• the muon acceptance, including the efficiency of the muon chambers, of the recon- 
struction and of the B/RMUON matching to central tracks, is known to about 7%. 
An independent dimuon sample was analysed to determine this uncertainty based on 
a method [39] which has been repeated here (±7%); 

• the error due to the uncertainty of the energy scale of the CAL was evaluated by 
varying the energy of the jets and the inelasticity i/jb in the MC by ±3% (±4%); 

• the efficiency of finding a track with 4 MVD hits was measured in the data and in 
the MC. The ratio of the measured efficiencies was applied as a correction to the 
acceptance. The uncertainty on this ratio was included in the systematic uncertainty 
(±3%); 

• the efficiency of the dijet trigger in the MC was corrected so that it reproduced the 
efficiency as measured in the data. The systematic uncertainty due to this correction 
was negligible; 

• the MC if" distribution was reweighted in order to account for the differences (see 
Fig. [3D between data and MC (< 1%); 

• the uncertainty on the size of the correction to the shape of the impact-parameter 
distribution for the MC samples described in Section H] was evaluated by varying the 
widths of the Gaussian and Breit-Wigner distributions used in the correction function 
by +20% and —10% of their nominal values. These variations are such that the global 
MC distribution still provides a good description of the data C^q%) ' 

• the uncertainty on the p™ 1 shape of the LF and charm background was evaluated by: 

— varying the correction applied to the LF background by ±20% of its nominal value 

(±2%); 

— varying the p^ 1 shape of the charm component by removing or doubling the cor- 
rection (±4%); 

• the contribution of flavour-excitation events in Pythia was varie ciby + 100%/- 50% 
and simultaneously the contribution of gg — > bb, qq — > bb events was varied by 
— 50%/ + 100% (±4%); the contribution of jg —> bb processes in Pythia was de- 
creased by 20% and all other processes were increased by +100% (±2%). 

The total systematic uncertainty was obtained by adding the above contributions in 
quadrature. A 2.6% overall normalisation uncertainty associated with the luminosity 
measurement was not included in the systematic uncertainty on the differential cross 
sections. 

3 The comparison between MC prediction and data for da/dxij* , see Section [3 is still satisfactory for 
these variations. 
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7 Results 



Figure [3] shows the distributions of the kinematic variables of the muon pt^ and rj^ as well 
as those for the jet associated with the muon p^T J and r]^ -3 . The fraction xZj* of the total 
hadronic E - p z carried by the two highest-p T is given by: 

jj _ E j= i, 2 (^ J ~Pz) , , 

X "~ E- Pz ■ {6) 

The distribution of x^Jj is also shown in Fig. [3l The data are compared in shape to 
the MC simulations in which the relative contributions of beauty, charm and LF were 
mixed according to the fractions measured in this analysis as described in Section HI The 
comparison shows that the main features of the dijet-plus-muon sample are reasonably 
well reproduced by this MC mixture. 

Total and differential visible cross sections have been measured for final states with at 
least one muon and two jets in the following kinematic region: 

• Q 2 < 1 GeV 2 and 0.2 < y < 0.8; 

• p 3 ^ 3 > 7, 6 GeV and \r] jl ' j2 \ < 2.5; the jets are defined as hadron- level jets using the kx 
algorithm. For the purposes of jet-finding, B hadrons are treated as stable particles; 

• p^> 2.5 GeV (p£ > 1.5 GeV for da/dp^) and -1.6 < rf < 1.3; 

• at least one muon is associated with a jet with p 3 T > 6 GeV. The muon is associated 
with the jet if it is the decay product of a B hadron contained in the jet, according 
to the hx algorithm. Muons coming from both direct (b — > //) and indirect (b — > 
c, c, J/ty, \&' — > /i) decays are considered to be part of the signal. 

The total visible cross section is 

a(ep -> ebbX -> ejjuX') = 38.6 ± 3.5(stat.) l^(syst.) pb. (4) 

This result is compared to the NLO QCD calculation described in Section [51 The predic- 
tion for the total visible cross section is 

a(ep -> ebbX -> ejjfiX') = 37.0^ 9 pb, (5) 
in excellent agreement with the data. 

4 xj? is the massive-jets analogue of the x° bs variable used for massless jets in other ZEUS publications 
[12]. 
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Figure H] and Table [TJ show the visible differential cross sections as a function of the muon 
transverse momentum and pseudorapidity. Also shown in Fig. H] and in Table [2] are the 
visible differential cross sections measured as a function of the transverse momentum of 
the jet associated with the muon p^~ 3 , and as a function of its pseudorapidity, rj^ 3 . The 
visible cross section as a function of pt^ is measured in the range pj> > 1.5 GeV, while the 
other cross sections are measured for pj, > 2.5 GeV. The NLO QCD predictions describe 
the data well and the Pythia MC also gives a good description of the shapes. 

The visible differential cross section as a function of rj^ is also compared with a previous 
ZEUS measurement [1], which used the p™ 1 method to extract the beauty fraction. The 
two measurements agree welfl. 

Figure [5J(a) and Table [3] show the visible dijet cross section as a function of xjj (Eq. [3]). 
The xif variable corresponds at LO to the fraction of the exchanged-photon momentum 
entering the hard scattering process. In photoproduction, events can be classified into two 
types of process in LO QCD. In direct processes, the photon couples as a point-like object 
in the hard scatter. In resolved processes, the photon acts as a source of incoming partons, 
one of which takes part in the hard scatter. The xij* variable provides a tool to measure 
the relative importance of direct processes, which gives a peak at x^ ~ 1, and of resolved 
processes, which are distributed over the whole x^ range. The dominant contribution to 
the visible cross section comes from the high-a^ peak but a low-x^ component is also 
apparent. The NLO QCD prediction describes the measured visible cross section well. 
Pythia also gives a good description of the shape of the distribution. 

Dijet angular correlations are particularly sensitive to higher-order effects and are there- 
fore suitable to test the limitations of fixed-order perturbative QCD calculations. The 
dijet variable measured, was reconstructed from the two highest-pr jets as: 

A0-" = -0 j2 | . (6) 

In the FMNR program, at LO the differential cross section as a function of A0-" is a 
delta function peaked at ir. At NLO, exclusive three-jet production populates the region 
|/T < < 7r, whilst smaller values of A0- 7 - 7 ' require additional radiation such as a fourth 
jet in the event. An NLO QCD calculation can produce values of Aft 3 ' < |-7r when the 
highest-pr jet is not in the accepted kinematic region. 

The visible differential cross section as a function of Aft 3 ' is shown in Fig. [5(b) and Table 
HI The NLO QCD predictions describe the data well. Visible cross sections as a function 

5 The measurement presented by HI [7] refers to a slightly different definition of the cross section and 
therefore cannot be compared to directly. However a qualitative comparison does not confirm their 
observation of an excess at low p!^. 
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of Afti have also been measured separately for direct-enriched (pjj > 0.75) and resolved- 
enriched {xjj < 0.75) samples (Fig. [5|c) and (d)). The cross sections are well described 
by the NLO QCD prediction for x{j > 0.75 and for xij> < 0.75. The Pythia MC gives 
an equally good description of the shape of the distributions. 

8 Conclusions 

Beauty production identified through semi-leptonic decay into muons has been measured 
with the ZEUS detector in the kinematic range defined as: Q 2 < 1 GeV 2 ; 0.2 < y < 0.8; 
p£' j2 > 7, 6 GeV; |r/ jlj ' 2 | < 2.5; f T > 2.5 GeV; -1.6 < rf < 1.3 with at least one muon 
being associated with a jet with pP T > 6 GeV. Lifetime information was combined with 
the muon p™ 1 method to extract the fraction of beauty events in the data sample. Unlike 
the previous analysis, which used the p™ 1 method alone, it was not necessary to constrain 
the charm contribution to the charm cross section obtained from other measurements. 
The extracted charm contribution is consistent with expectation. 

The total visible cross section was measured as well as visible differential cross sections 
as a function of the transverse momenta and pseudorapidities of the muon and of the 
jet associated with the muon. The rf cross section was compared to the previous mea- 
surement. This analysis confirms the previous result with similar statistical precision and 
different sources of systematic uncertainty. Also, it was possible to measure the cross 
section as a function of the muon transverse momentum to pt^. > 1.5 GeV, a lower pj, than 
in the previous muon-jet analysis. The p^, cross section agrees well with the NLO QCD 
prediction and does not confirm the excess observed by HI at low pl^,. 

All results were compared to the Pythia MC model and to an NLO QCD prediction. 
The NLO QCD prediction describes the data well. The Pythia MC model also provides 
a good description of the shape of the distributions. 

Beauty dijet angular-correlation cross sections were also measured. Separate measure- 
ments in direct-enriched and resolved-enriched regions were presented. The NLO QCD 
prediction describes the measured cross sections well. 
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range 


da/dpf^ ± stat. ± syst. 


^had 


(GeV) 


(pb/GeV) 




1.5, 2.5 


41.05±7.74 + H? 

— O.Ol 


0.87 


2.5, 4.0 


15.78 ± 1.96^-03 


0.93 


4.0, 6.0 


4.87 ± 1.03+°;^ 


0.98 


6.0, 10.0 


0.84 ± 0.27lg:n 


1.01 


r/ M range 


da/dr] fl ± stat. ± syst. 
(pb) 


Chad 


-1.60,-0.75 


3.86 ± 1.371^° 


0.83 


-0.75, 0.25 


16.81 ± 2.30+J;^ 


0.89 


0.25, 1.30 


19.70 ± 2.43+J;^ 


0.92 



Table 1: Differential muon cross section as a function ofpj, andr/^. For further 
details see text. The multiplicative hadronisation correction, Chad? applied to the 
NLO prediction is shown in the last column. 



Pt J range 
(GeV) 


da/dp^ 1 ± stat. ± syst. 
(pb/GeV) 


Chad 


6, 11 
11, 16 


4.74 ± 0.57ig;jjg 
1.78 ± 0.32+°;^ 


0.89 
0.89 


16, 30 


0.33 ± 0.10±g;gi 


0.92 


rf 1 '^ range 


da fdrf'i ± stat. ± syst. 
(pb) 


Chad 


-1.6,-0.6 


6.13 ±1.4lt^° 


0.77 


-0.6, 0.4 


13.89 ± 2.20±1;^? 


0.84 


0.4, 1.4 


16.42 ±2.29+^;™ 


0.99 



Table 2: Differential cross section for jets associated with a muon as a function 
°f Vt 3 &ndrf~ 3 . For further details see text. 
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xij range 


da/dxlj ± stat. ± syst. 
(pb) 


Chad 


0.000, 0.250 


11.85± 4.96 tlH 


0.69 


0.250, 0.375 


17.17± 7.89 t 7 A f 7 


0.78 


0.375, 0.500 


14.81± 7.56 tlf e 


0.86 


0.500, 0.750 


22.19± 4.48 tlH 


0.86 


0.750, 1.000 


106.63±12.63 t\lf 4 


0.92 



Table 3: Differential cross section as a function of x^ . For further details see 
text. 



range 


da/dAcf)" ± stat. ± syst. 
(Pb) 


Chad 


67r Sir 
12' 12 
8vr IOtt 
12' 12 

IOtt IItt 
12 ' 12 

117T 12tt 
12 ' 12 


2.26± 1.44 tUt 
7.35± 2.06 t\H 
24.70± 6.04 tUt 
92.91±11.10 tlUt 


0.80 
0.79 
0.86 
0.92 


xij > 0.75 


A<p^ range 


da/dA^i ± stat. ± syst. 
(Pb) 


Chad 


6-7T 107T 

12' 12 
IOtt 117T 

12 ' 12 
IItt 12tt 

12 ' 12 


1.62± 0.73 to™ 
15.27± 4.75 tHi 
65.69±10.66 tlii 


0.82 
0.87 
0.93 


xij < 0.75 


range 


da/dAcf)^ ± stat. ± syst. 
(Pb) 


Chad 


67r Sir 
12' 12 
8vr IOtt 
12' 12 

IOtt 117T 
12 ' 12 

IItt 12tt 
12 ' 12 


1.38±0.92 toil 
3.36±1.60 toil 
7.75±3.37 tlH 
18.84±4.17 tlil 


0.75 
0.76 
0.84 
0.84 



Table 4: Differential muon cross section as a function of A<fP 3 for all x^ and for 
x'ij > (<)0.75. For further details see text. 
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Figure 1: Distributions of (a) p™ 1 and (b) 5. The data are compared to a mixture 
of beauty (shaded histogram), charm (dotted line) and light flavour (dashed line) 
Pythia MC samples, combined with the fractions given by the two-dimensional 
Pt 1 -S fit. The total MC distribution is shown as the solid line. 
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Figure 2: Contours of 68% probability in the plane defined by the beauty and the 
charm fractions. The result of the x 2 fit to the two-dimensional p™ l -5 distribution 
(thick solid line) and the one-dimensional distributions in p™ 1 (thin solid line) and 
6 (dotted line) are also shown. The diagonal line shows the boundary of the physical 
region in which the fractions of b, c and LF are positive. 
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Figure 3: Distributions of (a) p^, (b) rf , (c) p^ J , (d) t/*"- 7 and (e) xij . The 
data are compared to a mixture of beauty (shaded histogram) , charm (dotted line) 
and light flavour ( dashed line ) P ythia MC predictions, combined according to the 
fractions given by the two-dimensional p^ cl -5 fit. The total MC distribution is 
shown as the solid line. The kinematic region is restricted to pt^ > 1.5 GeV (pj, > 
2.5 GeV) for (a)((b)-(e)). 
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Figure 4: Differential cross section as a function (a) pj,, (b) rj^, (c) p^ J and 
(d) Tf-i for Q 2 < 1 GeV 2 , 0.2 < y < 0.8, p£' j2 > 7,6 GeV, \r) juj2 \ < 2.5, and 
— 1.6 < ?/* < 1.3. For i/ie cross section, the kinematic region is defined as 
Pt > 1.5 GeV and as p 1 ^ > 2.5 GeV for all other cross sections. The filled circles 
show the results from this analysis and the open circles show the results from the 
previous ZEUS measurement. The inner error bars are statistical uncertainties 
while the external bars show the statistical and systematic uncertainties added in 
quadrature. The band represents the NLO QCD predictions with their uncertainties. 
The Pythia MC predictions are also shown (dashed line). 
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Figure 5: Differential cross sections as a function of (a) x 33 and (b) A<f) 3J of the 
jet-jet system and A<f> for (c) direct- and (d) resolved- enriched samples for Q 2 < 
1 GeV 2 , 0.2 < y < 0.8, p£' h > 7,6 GeV , rj juj2 < 2.5, f T > 2.5 GeV and -1.6 < 
rj^ < 1.3. The inner error bars are statistical uncertainties while the external bars 
show the statistical and systematic uncertainties added in quadrature. The band 
represents the NLO QCD predictions with their uncertainties. The Pythia MC 
predictions are also shown (dashed line). 
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